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We review the current status of the Brookhaven muon g - 2 experiment, and it’s 
effects on the SUSY parameter space when combined with dark matter relic density 
bounds, b ^ and Higgs mass constraints. If the 2>a deviation of g - 2 from the 
Standard Model value is correct, these data constrain the mSUGRA parameter 
space strongly, i.e. 300 GeV ~ mii 2 ~ §50 GeV, and mp ( at fixed tan/3, Aq) 
is tightly constrained (except at very large tan/3). Dark matter detection cross 
sections lie within the range accessible to future planned experiments. A non- 
universal gluino soft breaking mass however can greatly reduce the lower bound on 
mi /2 (arising from the b ^ S'y constraint) allowing for relatively light neutralinos, 
while non-universal Higgs H 2 mass can lead to new regions of allowed relic density 
where the detection cross sections can be increased by a factor of 10 or more. 


1 . Introduction 

While current Tevatron and LEP measurements have not greatly con¬ 
strained the SUSY particle spectra, there are a number of quantities, which 
if accurately measured and if accurate theoretical calculations existed, could 
greatly limit the SUSY parameter space of a given model, and thus allow 
significant predictions of what might be expected at the LHC and what 
might occur in the next round of dark matter detector experiments. We 
consider here the following quantities: the muon g - 2, the light Higgs 
mass rrih, the h —>■ sj branching ratio and the amount of dark matter. 
If accurately determined, these would greatly restrict the SUSY parame¬ 
ter space for a variety of models. We will first examine these within the 
framework of mSUGRA models (with R parity invariance) and then show 
some non-universal models which could moderate somewhat the mSUGRA 
constraints. 


1 


2 


2. The Muon g - 2 Anomaly: The Saga Continues! 

We review the current situation with the muon g - 2 magnetic moment 
anomaly. Recall that in 2001, the Brookhaven E821 experiment reported 
their high precision measurement of = {g — 2)/2. Based on the then 
best theoretical calculation of the Standard Model value ^, they reported a 
2.6cr deviation. Unfortunately, a sign error was subsequently found in the 
“scattering of light by light” (LbL) diagrams, which reduced the effect to 
only l.Gsigma. Since that time the following has happened: 

(i) New E821 data (the 2000 data) has been analyzed reducing the 
experimental error in by a factor of two^. The current world average is 
now 


a;fP = 11659203(8) x 10-^° (1) 

i. e. a measurement at the level of O.Tppm! One interesting feature of the 
Brookhaven measurements is the stability of their central value (with the 
error flag successively being reduced). 

(ii) New data has come from Novosibirsk (CMD-2) and Beijing (BES) on 
e^+e” ^ hadrons, and from ALEPH and CLEO on tau decay into two and 
four pions. These may be used to calculate the hadronic contributions to 
the SM prediction to get a more accurate determination of any deviation 
between theory and experiment that may exist. 

Two groups have now used this new data to reevaluate and we 
briefly review their results, can be divided into the following parts: 

= + + (2) 
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The QED and weak contributions to are well established, and the 
higher order (HO) hadronic contribution appears to be in good shape. With 
the corrected sign, the light by light (LbL) contribution has been evaluated 
by several groups ® with general agreement. We use here the value = 
[8.6 ± 3.5] X 10“^°. Current difficulties arise from the leading order in 
alpha (LO) hadronic contributions. This quantity can be calculated using 
a dispersion relation: 


a“ = 


^ 0 ) 




f Arni. 


ds^^R{s) 


( 4 ) 


where K{s) is the QED kernel and R{s) = CT(e“''e“ ^ hadrons)/cr(e+e“ ^ 
The integral is strongly weighted at low energy, i.e. about 90% 
comes from ^/s < 1.8GeV and about 75% from e+e“ ^ tttt (from the p). 
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Two procedures have been used to evaluate the dispersion integral. The 
hrst method uses the e+e" cross sections of CMD-2, BES and a large 
amount of earlier data to calculate R(s). However, it should be noted that 
the (T(e+e“ ^ hadrons) to be used in Eq.(4) is a “bare” cross section, 

1. e. the experimental data must be corrected for initial state radiation, 
photon vacuum polarization and electron vertex loop effects (and not car¬ 
rying this out correctly has led to some errors in past analyses.) Carry¬ 
ing out this analysis, Ref.^ finds a discrepancy between the experiment 
and the SM of Aa^ = 33.9(11.2) x 10“^° a 3 .Oct effect, while Ref.^ finds 

= 35.13(10.63) X 10“^° i.e a 3.3 ct effect. Thus the two analyses give 
consistent results. 

The second method makes use of the tan decay data of ALEPH and 
CLEO into 27r and 47r final states. Using CVC, the isovector form factor 
can be used to construct CT(e+e“ ^ 27r,47r) for s ~ 3GeV (which is most of 
the important region). However, in this case, one must include corrections 
due to the breaking of CVC. Major contributions to this come from the 
TT mass differences and the short distance radiative corrections (treated by 
chiral perturbation theory). Carrying out this analysis, Ref.^ finds from 
the tau data a discrepency of 16.7(10.7) x 10“^° i. e. only a 1.6 ct effect. 
Further, one cannot simply average the two proceedures of calulating ajj®. 
To point up the problem ^ reversed the proceedure, and using the e+e" data 
plus the CVC breaking corrections, predicted the tau branching ratios into 
27r and 47r final ststes and then compared this to the experimental ALEPH 
and CLEO data. They found that the prediction fails for the 7r“7r° state 
by 4.2 ct, and fails for the 27r“7r+7r° state by 3.5 ct. Thus the two approaches 
are statistically inconsistent with each other. 

At this point, there is no explanation for the disagreement between 
the two types of analysis. The discrepancy appears to be too large to be 
attributed to lack of understanding of the CVC breaking effects. Thus 
the question of whether the muon magnetic moment anomaly implies new 
physics is once again unclear. There will be more e“'"e“ data from CMD- 

2. BES and also KLOE and BABAR. The B-factories may also be able 
to measure the tau decays. Finally, the Brookhaven E821 experiment has 
3 X 10® /i“ events it is currently analyzing. (Results may be out by early 
next year.) It also has future plans for running the experiment further. 

3. mSUGRA Model 

In SUSY models, one generally has a contribution to in addition to 
the SM piece This arises from loops with chargino and sneutrino or 
neutralino and smuon intermediate states, with a magnetic field attached 
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Figure 1. mSUGRA allowed regions for tan/3 = 10, Aq = 0, /r > 0. The dot-dash lines 
are the bounds for a Linear Collider(LC) Xi ~ X 2 signal (for ^/s = 500 GeV and 800 
GeV), and the curved lines are for the stau pair production signal. 

to any charged particle. In general these contributions are not small. If the 
analysis leading to a 3cr effect is valid, and if one attributes the deviation 
to the SUSY correction, one gets a significant constraint on the SUSY 
parameter space, with the SUSY spectrum required to be relatively low 
and easily within the reach of the LHC. On the other hand, if the 1.6 sigma 
analysis turns out to be correct implying only a small SUSY contribution 
is possible, then squark and gluino mass spectrum would be pushed to the 
TeV domain. Thus the resolution of the current ambiguity in is very 
important. In this section, we analyze these matters within the framework 
of the mSUGRA model with R parity invariance. 

mSUGRA is the simplest SUGRA model in that it depends only on four 
parameters and one sign: mo (the scalar soft breaking mass at Mg), mi /2 
(the gaugino mass at Mq), Aq (the cubic soft breaking parameter at Mq), 
tan/3 =< H 2 > / < Hi > (at the electroweak scale), and the sign of /i 
(the Higgs mixing parameter in the superpotential:/iiLiiL 2 ). We examine 
here the parameter range mo > 0, mi /2 < ITeV (which corresponds to the 
gluino mass bound accessible to the LHC: nig < 2.5GeV), 2 < tan/3 < 55, 
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Figure 2. Same as Fig.l for tan/3 = 40, Aq = 0, /i > O.The dashed line give contours 
of Bs —> /i/i branching ratios at the Tevatron. The vertical dotted lines are Higgs mass 
contours. The short vertical lines are cr^Q_^ = 3 x 10“^pb (lower) and 1 x 10“®pb 
(upper). 


and l^ol < 4 toi/ 2 - One starts the analysis at Mq and uses the renormal¬ 
ization group equations (RGE) to go down to the electroweak scale. Thus 
all SUSY masses and cross sections are determined in terms of these four 
parameters and one sign. The details of carrying out this analysis including 

all coannihilation effects in the relic density analysis can be found in e. g. 
8 

The mSUGRA model allows one to calculate an array of quantities that 
can be measured and hence can be used to restrict the parameter space. 
The include the neutralino relic density, the 5 —> s -I- gamma branching 
ratio and the Higgs mass bound. For the relic density we take a 2cr range 
around the current central value 

0.07 < < 0.21 (5) 

For the 5 —> s -I- 7 decay (which has both systematic and theoretical uncer¬ 
tainties) we take a relatively broad range around the CLFO central value 
and use the LEP bound for mh 

1.8 X 10“^ < B{b sj) < 4.5 X 10“^; > IMGeV 


( 6 ) 
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In addition there is the LEP bound on the chargino mass of > 

103GeVi2. 



mj/2[GeV] 

Figure 3. Same as Fig. 2 for tan (3 = 55, Aq = 0, /i. > 0 and mt = 175GeV, mi) = 4.25. 

In spite of the large errors in the experimental data that still are present, 
Eqs.(5,6) put significant constraints on the SUSY parameter space in that 
they require that TOi /2 (300 — 400) GeV across the full parameter space. 
This is illustrated in Figs. (1-3) for tan/3 = 10, 40, 55. If we assume the 
3cr deviation for the muon is valid, then the 2a bound from the central 
value is 11 x 10“^°. For low tan/3, e.g. tan/3 = 10 (Fig. 1), the combined 
constraints then leave very little available parameter space. For higher 
tan/3, e.g. tan/3 = 40 (Fig.2), more allowed region exits, but the xt mass 
constraint combined with the constraint eliminates all the “focus point” 
region of large toq and low toi/ 2 . For very high tan/3, e.g. tan/3 = 55 
(Fig.3) a new region opens up producing a “bulge” at low mi /2 due to 
rapid s-channel annihilation of the A and H Higgs bosons (as their mass is 
significantly reduced). Throughout the entire range, the Xi ~ proton cross 
sections for direct detection of Milky Way dark matter lie in the range of 
(9 X 10“® — 5 X 10“^°)pb, a range that is expected to be accessible to future 
large scale experiments. 
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Figure 4. Same as Fig. 2 for tan/3 = 40, Aq = M > 0 non-universal soft 

breaking of the Higgs masses with ^2 = 1. The lower shaded band is the usual allowed 
stau-neutralino coannihilation band, and the upper band is the new region arising from 
the non-universal Higgs masses due to increased annihilation through the Z-channel. 



4. Non-Universal Models 

We consider here two types of non-universal soft breaking at the GUT scale: 
for the gaugino masses and for the Higgs masses. The first is of interest 
in that it can soften significantly the lower bound produced by the 5 —> sy 
and rufi constraints. Thus if one assumes at Mqut that the gluino mass is 
7711/2(1 + S 3 ) (where 7711/2 is the universal mass), then for i 53 = 1 , one finds 
e.g. for tan/3 = 50,that the lower bound on ?77i/2 is reduced to 185 GeV 
corresponding to a neutralino mass of 75 GeV. (The constraint from rrih is 
also softened, though the constraint becomes somewhat stronger). 

If the H 2 mass is increased at Mq, new effects also occur. Thus writing 
^772 “ ^^ 0(1 + '^ 2)1 where mo is the universal scalar mass, then a new 
annihilation channel in the relic density analysis can open up for low 7711/2 
and high mo due to rapid annihilation through the s-channel Z pole. This 
is shown in Fig. 4 for tan /3 = 40. In this new region, the neutralino-proton 
cross section also increases by a factor of 10 or more, allowing cross sections 
in the range (10“® — 10“^)pb. 













5. Conclusions 

The muon magnetic moment anomaly can produce strong constraints on the 
allowed region of SUSY parameter space. Thus if the 3a deviation analysis 
is correct, then a 2 a bound from the central value gives an upper bound 
of TOi /2 ~ 850GeV when combined with the relic density constraint. The 
b —>■ s"f and mh constraints in the mSUGRA model, produce a lower bound 
of TOi /2 > (300 — 400) GeV, thus bounding the parameter space in a region 
easily accessible to the LHC. Dark matter detection cross sections then lie in 
a region accesssible to future dark mater experiments.However, should the 

anomaly be smaller than ~ 10 x 10“^° the squark and gluino spectrum 
will be pushed into the TeV domain. Non-universal models can modify these 
constraints. Thus an increase in the gluino mass at Mq can signficantly 
decrease the the lower bound on TOi /2 (and hence on the neutralino mass), 
while an increase of the H 2 mass at Mq can give rise to new allowed regions 
of relic density with detection cross sections increased by a factor of 10 or 
more. 

This work was supported in part by the National Science Foundation 
Grant PHY - 0101015. 
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